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ABSTRACT 

Deep surveys in many wavebands have shown that the rate at which stars were forming 
was at least a factor of 10 higher at redshifts > 1 than today. Heavy elements ('metals') 
are produced by stars, and the star formation history deduced by these surveys implies 
that a significant fraction of all metals in the universe today should already exist at 
z ~ 2 — 3. However, only 10% of the total metals expected to exist at this redshift 
have so far been accounted for (in Damped Lyman Alpha absorbers and the Lyman 
forest). In this paper, we use the results of submillimetre surveys of the local and high 
redshift universe to show that there was much more dust in galaxies in the past. We 
find that a large proportion of the missing metals are traced by this dust, bringing 
the metals implied from the star formation history and observations into agreement. 
We also show that the observed distribution of dust masses at high redshift can be 
reproduced remarkably well by a simple model for the evolution of dust in spheroids, 
suggesting that the descendants of the dusty galaxies found in deep submm surveys 
are the relatively dust-free spiral bulges and ellipticals in the universe today. 

Key words: galaxies: evolution - galaxiesiISM ~ galaxies: luminosity function, mass 
function - submillimetre ~ (ISM):dust, extinction 
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1 INTRODUCTION 

The star formation history (SFH) of the universe deduced 
from optical/UV surveys (Lilly et al. 1996; Madau et al. 
1996, 1998) has recently been revised to account for dust 
extinction (Steidel et al. 1999), largely motivated by the 
detection of a population of dusty star-forming objects in 
deep submm surveys with SCUBA (Hughes et al. 1998; Bales 
et al. 1999; Smail, Ivison & Blain 1997; Scott et al. 2002). 
The current best estimate of the SFH is that the global 
star-formation rate in the universe was fairly constant over 
the redshift interval 1 < 2 < 4, falling by a factor > 10 to 
the present epoch, see Fig. 1. (Steidel et al. 1999; Ivison et 
al. 2002). A simple integration of this history implies that 
~ 25% of all stars and metals today should have formed by 
a redshift of ~ 2.5.* However, the metal content of objects 
believed to contain most of the baryonic matter at high-z 
- principally the Lyman alpha forest and damped Lyman 
alpha absorption systems (DLAs)- is only 2.5% of the total 
amount of metals in the universe today (Pettini 1999; Pagel 



In a flat cosmology. 



2002). Thus 90% of metals expected to exist at 2 = 2 — 3 
have not yet been seen or recognised - the 'missing metals'. 

The submm surveys revealed a population of extremely 
dusty and obscured objects with very high bolometric lumi- 
nosities, probably residing at high redshifts, z > 2 (Dunlop 
2001a; Ivison et al. 2002). If the dust is heated by stars, 
the star formation rates implied by the submm fluxes are 
~ 100 — lOOOM0yr~^, high enough to produce an entire 
galaxy worth of stars in ~ 1 Gyr. Since the properties of low- 
redshift ellipticals suggest that most of their stars formed 
quickly at high redshift (Bower et al. 1992; Renzini 1997; 
Jimenez et al. 1999), several authors have argued that the 
SCUBA sources may be elliptical galaxies in the process 
of formation (Lilly et al. 1999; Bales et al. 2000; Dunlop 
2001b; Granato et al. 2001). It has been claimed (Cen & 
Ostriker 1999) that the regions of greatest over-density in 
the universe at any epoch are the most metal enriched, and 
that these regions will have reached a given metallicity much 
earlier than lower density environments. This argument sug- 
gests that SCUBA galaxies should be a large reservoir for 
metals at high redshift, since if they are proto-ellipticals they 
should have formed from the highest density peaks at high 
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Figure 1. The current best estimate of the star formation history 
of the universe, taken from Ivison et ah 2002. The points repre- 
sent extinction corrected optical/UV, radio and submm estimates 
while the Unes are the latest submm estimates based on the 8mjy 
survey (Scott et al. 2002) using different redshift estimators. The 
arrow represents the extrapolation of the 8 mjy population to 
1 mJy (see Ivison et al. 2002). All estimates are in reasonable 
agreement, leading to a star formation history which is relatively 
constant between z = 1 — A, declining by a around a factor of 10 
from z = 1 to today. 



redshift when there has been less time for pollution of the 
low density inter-galactic medium (IGM). 

In this paper, we wish to investigate the connection be- 
tween the 'missing metals', the formation of the spheroids 
and the dusty sources seen by SCUBA. However, instead 
of using the bolometric luminosities of the submm sources 
to infer their physical properties, we interpret the results 
of the deep submm surveys and our submm survey of the 
nearby universe in a completely different way, by carrying 
out the first investigation into the cosmic evolution of the 
dust mass in galaxies. This provides a new way of looking 
at the chemical evolution of the universe, as there is evi- 
dence that the fraction of ISM metals which are bound up 
in dust grains is ~ 40 percent, a value which seems remark- 
ably constant from galaxy to galaxy and which also appears 
to apply at higher redshifts (Issa, MaClaren & Wolfendale 
1990; Edmunds 2001; Pei, FaU & Hauser 1999; James et al. 
2002). We will determine the dust content of the low and 
high-z universe by constructing dust mass functions, i.e. the 
space densities of galaxies as a function of dust mass. At 
low-z, we will use the results of the SCUBA Local Universe 
Galaxy Survey (Dunne et al. 2000), the first attempt at an 
unbiased submm survey of the local universe. For the high 
redshift dust mass function, we will use submm data from 
the deep SCUBA submm surveys. In Section 4 we will use 
these dust mass functions, together with a simple chemical 
evolution model, to estimate the co-moving density of met- 
als and baryons associated with the ISM of galaxies at low 
and high redshift. Finally, we will produce a simple predic- 
tion for the dust mass function of the spheroids during the 
epoch in which they had their maximum dust masses, and 



compare with our observational estimate of the dust mass 
function of the SCUBA sources. 



2 DUST MASSES FROM SUBMILLIMETRE 
FLUXES 

Dust mass can be calculated from the submm flux according 
to 



(1-l-z) D'' So 



(1) 



Kd(j^em) B(:/em, (Td)) 

where D is the distance for the assumed cosmology and 
Sobs is the observed flux. The dust mass opacity coefficient 
^d(iycm) normalises the submm emission for a given dust mass 
and is currently measured only at FIR wavelengths and then 
extrapolated to the submm with a A"'' dependency. We 
have taken the average value from the following sources of 
«:d(i25fim) ~ 2.64 ± 0.29 m'^ kg~^, where the ±0.29 is the er- 
ror on the mean value (Hildebrand 1983; Draine & Lee 1984; 
Kriigel, Steppe & Chini 1990; Sodroski et al. 1997; Bianchi, 
Davies & Alton 1999). We use (3 = 2 (for which there is com- 
pelling evidence from our local survey and elsewhere, Dunne 
& Eales 2001 and references within) to extrapolate Kd(i25Aim) 
to the rest wavelength of the submm emission (~ 240/im for 
an 850^m source at z ~ 2.5). This agrees with the recent es- 
timation of Kti(850Mm) = 0.07±0.02 w? kg^^ using a different 
technique (James et al. 2002). 

We also require the value of {Td), which is a mass- 
weighted dust temperature. Submillimetre and ISO observa- 
tions of many nearby galaxies have shown that, when enough 
data is present, the SED is best described by dust at more 
than one temperature (Haas et al. 1998,2000; Frayer et al. 
1999; Braine et al. 1997; Dumke et al. 1997, Calzetti et al. 
2000; Guelin et al. 1995; Sievers et al. 1994; Neininger et 
al. 1996; Aflon et al. 1998, 2001; Papadopoulos & Seaquist 
1999; Reach et al. 1995; Trewhella et al. 2000). For our local 
SCUBA survey of bright IRAS galaxies, we found that in 
addition to the typical TRAS-detected' dust at 35 — 45 K, 
a colder dust component at 17 — 25 K was required in all 
cases (Dunne & Eales 2001), even for highly luminous galax- 
ies such as Arp 220. Klaas et al. (2001) also find that the 
SEDs of 41 local ULIRGs are best represented by a combi- 
nation of warm and cold dust. Following the notation of the 
two-component SED model from Dunne & Eales (2001), the 
value which should enter Eqn. 1 can be approximated by 

2(7ui Xc 



(2) 



where T^jT^ is the characteristic temperature of the 
cold/ warm dust, and Nc/N^ is the ratio of mass in the 
cold/warm components. For the local survey, the value of 
{Td) was 21.3 K, with a standard error in the mean for 
the population of 0.5 K. This mass-weighted temperature 
should not be confused with the dust temperature tradition- 
ally used to estimate bolometric luminosity, which should be 
the temperature defining the peak in the spectrum. Since we 
are not integrating the SED and have observations at wave- 
lengths longer than the peak, our estimation of the dust 
mass from a single submm fiux is not as sensitive to tem- 
perature as an estimate of the bolometric luminosity would 
be. In local galaxies, including ULIRGs, most dust is con- 
tained in the cold ~ 20 K component, which means we do 
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not need to worry about the location of the peak in the SED; 
the temperature of the warm dust has a neghgible effect on 
the dust mass estimate, providing Nc/Nw > 1. For the local 
objects, we use (Td) for each galaxy in our survey (Dunne 
& Eales 200 i). 

The biggest uncertainty in our analysis is our lack of 
knowledge of (Td) in the high redshift SCUBA galaxies. As 
discussed above, there is now a lot of evidence that most 
of the dust in low redshift FIR sources, including ULIRGs, 
is at a relatively cold temperature. It will not be possible 
to perform a similar analysis of the temperature of dust in 
high redshift galaxies because of the lack of instruments op- 
erating at wavelengths long enough to sample the rest-frame 
emission from the cold dust (i.e. 850/im in the rest-frame). 
It is difficult to determine the existence of a colder dust com- 
ponent without many fluxes at wavelengths at and beyond 
the peak in the SED. It is plausible that the mean tem- 
perature of dust in galaxies was higher in the past, as star- 
formation rates and possibly interstellar radiation fields were 
higher. However, it is also possible that the mass-weighted 
dust temperatures in high-z SCUBA galaxies might be as 
low, or even lower than in low-z ULIRGs. This is because 
dust equilibrium temperatures are very sensitive to the rela- 
tive geometry of dust and heating sources, as well as optical 
depth effects. If the clouds in which the stars form are very 
dense, it is likely that few UV/optical photons escape (we 
know this is the case for the high redshift SCUBA galax- 
ies as a whole), and so dust outside these regions will see 
a much weaker radiation field and hence have a lower tem- 
perature. There is also evidence (Prayer 2001; Ivison et al. 
2001,2002; Isaak et al. 2002; Lutz et al. 2001, Eales ot al. 
2000) that the physical sizes of the star forming regions in 
some submm luminous objects at high redshift are larger 
than in local ULIRGs. Local ULIRGs have very compact re- 
gions of star formation (< 1 kpc) while the radio, mm and 
CO sizes of the high-z submm sources have been estimated 
to be as large as 10-30 kpc in some cases. This could be 
due to the high-z objects having much greater gas fractions 
than local ULIRGs. Mihos & Hernquist (1996) have argued 
that it is the presence of largo stellar bulges in local ULIRGs 
which provide enough stability to allow the gas to reach very 
high central concentrations before forming stars. Galaxies 
at high redshift may have greater gas fractions and small 
or non-existent stellar bulges compared to galaxies today. 
This means that the star formation induced by an interac- 
tion or perturbation at high redshift could be much more 
widespread, as the gas will collapse to form stars before it 
inflows to the centre of the remnant (Mihos 1999; Tissera et 
al. 2002). A given star formation rate which is spread over a 
larger volume will result in a lower intensity radiation field, 
possibly leading to lower dust temperatures. Two sources at 
« = 0.5 - 1 from the 170/jm FIRBACK survey have also 
been followed up with SCUBA and appear to have much 
colder dust for their luminosity than local objects selected 
by IRAS (Chapman et al. 2002). Such sources may be more 
representative of objects selected in deep submm surveys, 
as selection at longer wavelengths is less biased toward ob- 
jects dominated by warmer dust. The models of Pei, Pall & 
Hausor (1999) however, predict that there should have been 
a modest amount of positive evolution in the dust equilib- 
rium temperature between z = and x = 2 — 3. In the 
absence of any direct evidence on the amount of cold dust 



in the high-z SCUBA galaxies, we have assumed the equilib- 
rium temperature predicted for dust sources at ~ 2 — 3 by 
Pei, Pall & Hauser (1999) of 25 K. If temperatures are more 
akin to local galaxies (20 K) then the high-z dust masses 
we present will increase by a factor ~ 2, and more if the 
equilibrium temperatures at high redshift are actually lower 
than in local ULIRGs. 



3 DUST MASS FUNCTIONS AT LOW AND 
HIGH REDSHIFT 

We estimate the space density of galaxies as a function of 
dust mass (the dust mass function - DMF) using standard 
accessible volume techniques, such that 

$(M)AM = E,-i- (3) 

where Vi is the accessible volume in the parent sample of 
the ith source and the sum is over all objects in the range 
M — > M -h AM. For the local galaxies Vi is their accessible 
volume at 60/itm, which was the wavelength of selection of 
their parent sample, the IRAS BGS (Soifer et al. 1989). The 
local dust mass function (DMF) we will use here is described 
in detail in Dunne & Eales (2001). 

For the high redshift dust sources, we construct the 
DMF using submm data from the CUDSS blank field submm 
survey (Eales et al. 2000; Webb et al. 2002), and also the 
submm number counts at higher and lower flux densities 
(Small et al. 2002; Scott et al. 2002; Borys et al. 2002). In 
order to allow for star forming objects which may be too 
submm faint to appear in the blank field surveys, we use 
the statistical submm analysis of starburst galaxies in the 
HDP by Peacock et al. (2000). The characteristic shape of 
the thermal dust spectrum means that in a flat cosmology, 
an object of a given dust mass will be observed to have vir- 
tually the same 850/im flux-density over a wide range of red- 
shift {1 < z < 10). Therefore the redshift we assume for the 
sources has very little effect on the derived dust masses, pro- 
viding that the objects are at a > 1 (any objects known to 
have z < 1 were excluded from this analysis). Spectroscopic 
redshifts (Eales et al. 2000) and estimates using the Dunne, 
Clements & Eales (2000) version of the radio-submm red- 
shift indicator (Caxilli & Yun 1999) were used where avail- 
able, otherwise sources were placed at z = 2.5 for the pur- 
poses of calculating Md. Figure 2 shows that due to the 
flat flux-density redshift relation, a source visible at z = 1 
will still be visible (in a flux limited sample) at a = 5 — 10 
(depending on the exact shape of the spectrum); therefore 
Vi for the submm sources was simply taken to be the vol- 
ume in the redshift interval Zmax — Zmin- We fix Zmin = 1 as 
there is dramatic evolution in the submm population below 
this redshift, very few submm sources being found at « < 1 
(Dunlop 2001a), while we take 2max = 5. This is the red- 
shift at which the 850/xm flux density of a source with an 
SED similar to Arp 220 would start to fall off (Fig. 2). The 
effects of using ^max = 4 or ^max = 6 are a modest ±20% 
change in the space density. We have treated the statistical 
detection of submm flux in the HDP in a similar manner 
to the blank field sources, using z = 3, the average of the 
photometric redshifts, to estimate the dust mass. The pho- 
tometric redshifts show that the galaxies responsible for the 
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Figure 2. The 850/tm flux density for a source with a given 
dust mass in an Om = 0.3 , Ha = 0.7 cosmology. Various SEDs 
are shown: SLUGS single temperature is the average single tem- 
perature fit to the local survey data (/3 = 1.3, T^j = 36 K ) 
Dunne et al. (2000); average SLUGS 2-component fit provides 
the mass-weighted temperature of 21.3 K, [Tu, = 42 K, Tc = 21 
K, Nc/Nw = 10]; warmer SED with (T^) = 25 K as assumed for 
the high-z sources, [T™ = 45 K, Tc = 25 K, N^/N^ = 13], similar 
to that fitted to ULIRG IRAS 1211+03 (Dunne & Eales 2001); 
cold SED for more quiescent galaxies like NGC 891 (Alton et al 
1998), [r„ = 30 K, Tc = 15 K, Nc/N^ = 100]; the 2-component 
SED for Arp 220 (Dunne & Eales 2001) [T™ = 48 K, Tc = 18 K, 
Nc/Nw = 42], All 2-component SEDs have /3 = 2. The fluxes for 
the 2-component SEDs drop off after z = 4 — 5, most dramatically 
for those with more cold dust. 

submm emission are found between z = 1 — 6 and this red- 
shift range was used to calculate the space density. There 
were also 5 optical sources in the HDF which are plausible 
> 3 cr detections in the submm, these are treated separately 
and we have used their photometric redshifts to estimate the 
dust mass. 

This high-z dust mass function is shown in Fig. 3 
for two cosmologies, ilm = 1 (solid) and a A model with 
Qa = 0.7, flm = 0.3 (blue dashed), along with the DMF 
from our local survey. A value of Ho = 75kms~^ Mpc~^ is 
used throughout. It is immediately obvious that there has 
been substantial evolution in the dust mass function, par- 
ticularly at the highest dust masses. 



4 THE METAL CONTENT OF THE 
UNIVERSE 

There is growing concern that the metal content of the high- 
z universe is not observed to be as large as the inferred star 
formation would suggest (Pettini 1999; Pagel 2002; Blain et 
al. 1999a,b). Accounting for the metals observed in DLAs 
and the Lyman forest leaves a shortfall of ~ 90% in the cos- 



mic metal budget at 2; = 2 — 3. Metals are thought to be 
produced by supernovae and in the winds of cool evolved 
stars. They are then injected into the ISM with a fraction 
of them (77) bound into, or subsequently accreting onto dust 
grains. From here they can be expelled out of the galaxy 
into the inter-galactic medium (IGM) by starburst winds, 
or become part of the next generation of stars. The fraction 
of metals found in these different locations today is sum- 
marised in Fig. 4a. 

These fractions do not remain constant, of course, as 
the fraction of total metals in stars and the IGM will de- 
crease with redshift. Today the fraction of metals in the 
ISM of galaxies is negligible, but how different could this 
have been in the past, when galaxies were more gas rich? 
It is clear that more metals should be found in the most 
evolved objects at a given redshift, although their distribu- 
tion between stars and ISM will depend on how far the ob- 
ject has evolved. This has been presented by Cen & Ostriker 
(1999) as a correlation between over-density and metallic- 
ity, the fastest evolving objects being found preferentially in 
the most over-dense regions. The monolithic collapse models 
predict massive episodes of star formation at high redshift, 
and this should be accompanied by rapid enrichment of the 
locality. Observations of DLAs and the IGM at high redshift 
have failed to find evidence for this level of enrichment, al- 
though high-z QSOs are often found to have near solar or 
super-solar metallicities (Warner et al. 2002). This makes 
sense if QSOs represent the accretion phase of the massive 
black hole, which is presumably forming co-evally with the 
stellar bulge, as required to produce the BH-bulge mass cor- 
relation in the local universe (Magorrian et al. 1998). The 
IGM is a low-density environment and so not expected to 
be very enriched at high-z. It has been argued that the 
Damped Lyman alpha absorbers represent the co-moving 
mass density of gas in galaxies at high redshift and that the 
column density-weighted metallicity of these systems traces 
the global metallicity evolution of the universe^ (Pei & Fall 
1995). However, there are potential biases in the selection of 
DLAs as they are only observed along sight lines to optically 
bright QSOs and thus they are less likely to be observed if 
they contain dense gas and dust. The cross-sectional area 
available at the edges of galaxies is also greater than the 
denser, more evolved inner regions and therefore a sight-line 
is more likely to intercept a less dense and less metal en- 
riched part of a galaxy (Mathlin et al. 2001). Thus current 
samples of DLAs may underestimate the metal content of 
the neutral gas at any redshift. The obvious place to look 
for the most enriched systems at high-z is to look at the 
most obscured and actively star forming objects - submm 
sources, especially as they show evidence for being located 
in over-dense regions (Ivison et al. 2000). 

One simple thing we can do with the high-z DMF is to 
integrate the total amount of dust in galaxies at low and 
high redshift. The co-moving density of dust is given by 
Pd = r(2 -f a)M^(l)* , where and (j)* are the Schechter 
parameters fitted to the dust mass functions in Fig. 3. pd 
can be converted into pz, the co- moving density of met- 
als, by dividing by 7; - the fraction of ISM metals which 

t provided that no class of galaxy is excluded from the DLA 
samples 



© 0000 RAS, MNRAS 000, 000-000 



A census of metals at high and low redshift 5 




dust mass (M^) 

Figure 3. Dust mass functions at iiigli rcdsiiift (diamonds and triangles) and locally (circles) along with best fitting Schechter functions. 
The observations at high redshift match the predictions for spheroids (shaded region) when their dust masses were at a maximum. High-z 
points estimated using data from: CUDSS blank field SCUBA survey (Bales et al. 2000; Webb et al. 2002) — points 4,5,6; submm source 
counts at higher and lower fluxes (Small et al. 2002; Scott et al. 2002; Borys et al. 2002) - points 2,7; submm properties of optical 
starbursts in the Hubble Deep Field (HDF) Peacock et al. (2000) — points 1,3. Errors are simple Poisson statistics. Blue dashed line 
and open triangles are for a Qrn = 0.3, = 0.7 cosmology, solid line and points for = 1. Shaded region indicates the maximal dust 
evolution for local spheroids assuming (Af/L)* = 4.5 — 8.3 and a closed box chemical evolution model. The space densities for the high-z 
DMF change by ~ ±20% if the volume is calculated using 2max = 4, 6 rather than 5. Schechter parameters for the high-z DMF and 
Om = 1 are a = -1.22, M* = 2.7 x 10** Mq, 4>* = 2.4 x 10"^ Mpc-^ (M* increases to 5.4 X 10® Mq for (r<j> = 20 K) , and for the A 
cosmology are a = -1.08, M* = 4.7 x 10* Mq, <j>* = 8.9 x 10"** Mpc'^. The Schechter parameters for the local DMF are a = -1.85, 
M* = 9.3 X 10'' Mq, (f>* = 2.1 X 10-''Mpc-3. 



are in dust. We have used rj = 0.4 and there is some evi- 
dence that this fraction is fairly constant both locally and 
at higher redshifts (Pei, Fall & Hauser 1999; Issa, MacLaren 
& Wolfendale 1990; Edmunds 2001; James et al. 2002). Ta- 
ble 1 lists the ISM metals at low and high-z inferred from the 
dust, along with an inventory of metals and baryons in other 
locations. Integrating the DMF is currently quite uncertain 
at low redshift due to the very steep faint end slope fitted 
to the local The minimum and maximum values in 

Table 1 are obtained by integrating the local DMF down to 
the lowest observed mass, and to zero mass respectively. To 
facilitate comparisons with work from the literature, we now 
restrict ourselves to the i},n ~ 1 cosmology, noting that this 
does not affect any of our conclusions. 

The baryon values are taken from the literature 
(Fukugita, Hogan & Peebles 1998 - FHP) except for that 
associated with the ISM at high-z. For this we have used 

t This may simply be because we have not yet probed to low 
enough dust masses to properly sample the turn over 



the metals in dust, plus a simple 'closed-box' chemical evo- 
lution model to predict the baryon content of the ISM at 
high redshift. This chemical evolution model (Fig. 5), and 
others which incorporate inflow and outflow of gas, predict 
that the ratio of dust mass to baryonic mass should change 
quite slowly over a large part of a galaxy's evolution, and 
should peak when roughly half the gas has been converted 
into stars (Eales & Edmunds 1996, Edmunds & Eales 1998). 
We have used the maximum ratio of dust mass to baryonic 
mass (we are assuming that the submm sources are observed 
at their maximum dust mass) given by the model, and the 
integral of the DMF, to estimate the total amount of baryons 
associated with the objects in the DMF. Since the maximum 
occurs at a gas fraction of ~ 0.5, the total baryonic mass is 
split equally between stars and gas (ISM). The assumption 
that the submm sources are observed close to their peak 
dust mass is probably valid for the bright submm selected 
sources, but cannot be presumed to hold for the faintest 
lensed sources and the statistical detection of optically se- 
lected sources in the HDF. This means that the baryon den- 
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Metals at z=0 Metals at z=2.5 

1 ^nvin^i 




stars ISM DLA IGM IGM+ stars ISM DLA IGM IGM + 



(c) (d) 

Figure 4. The location of metals and baryons at low and high redshift. The categories are taken from Table 1. DLAs are given as a 
separate category because it is not clear whether they are produced by interstellar gas or intergalactic gas, or a mixture of both. The 
darker shaded IGM represents that currently observed - the hot cluster gas at low redshift and the Lyman forest clouds at high redshift. 
The lighter IGM+ represents the range of values of extra IGM which have been inferred for the group/field gas at low redshift and 
from CDM predictions for the total Lyman forest gas at high redshift (FHP). The unshaded regions for ISM and IGM+ in (a) and (c) 
represent the maximum which could be present within the uncertainties described in Table 1. For low-z ISM, this is when the integral of 
the DMF is extrapolated to zero mass. The additional IGM+ gas would appear to make up most of the baryonic content at all redshifts 
(without it, there is a serious shortfall in accounting for the baryon content expected from CMB analysis and nucleosynthesis), but is 
not an important location for metals at high redshift. Most metals at high-z are evidently in the ISM, while at low redshift they are in 
stars and, possibly, the low density IGM+. 
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Table 1. Baryons and Metals in the local and high redshift universe 



Location 



Baryons (M© Mpc ^) 
= 2 = 2.5 



Metals (MoMpc-^) 
z = 2 = 2.5 



stars 



ISM 



DLAs 



5 X 10* h76 ^ 

9.2 X 10'' h75 ^ 
H1 + H2 



IGM observed 3.7 x 10* h°-^ ^ 
cluster gas 

IGM+ 1 - 4.4 X 10^ h75 ' 

group/field gas 



1.1 X 10* h75 

/SFH 

1.3 - 2.7 X 10* h75 
DMF & model 

2.2 X 10* h75 a 

ISM / IGM? 

1.4 X 10^ h75^ 
Lya forest 

4.2 X 10^ h7B a 
extra inferred 
by CDM 



9.6 X 10^ h75 

Z = Zq'' 
0.8 - 3.0 X 10^ h75 



2.3 X 10® hOgS 
Z = 0.33^0 <^ 

1.9 - 8.3 X 10« h75 
Z = 0.1 Z0 



» 4.4 X 10^ h75 

sec text 

1.9-3.8 X 10** h75 
(Td) = 25 - 20 K 

2.9 X 10^ h75 

Z = 0.07Zq 

8.2 X 10* h75 
Z = 0.003Zq >> 

2.4 X 10^ h75 
Z = O.OO3Z0 



Total observed 
Inferred 



l-5.4xl0^h75 1.9 - 6.2 X 10** h76 



6.1 X 109< 



6.1 X 10^' 



1.2 - 2.1 X lO'' h75 2.7 - 4.9 x 10^ h7B 
1 - 2 X lO'' 2.5 - 5 X 106 



Table adapted from Pettini (1999) and Pagel (2002) with h76 = Ho/75kms~'^ Mpc~^ and Qm = 1 for consistency with previous work. 
We assume Zq = 0.0189. The metal content has been estimated from either the baryon content and the metallicity, or from the dust 
masses. The 'IGM observed' represents the directly observed gas (the hot cluster gas at low-z, and the Lya clouds at high-z). The 
'IGM+' is the additional IGM gas which has been inferred but not directly observed (the group/field gas locally, and predictions for 
the total IGM from CDM simulations at higli-z). The 'total observed' row has a minimum when counting only 'IGM observed' and a 
maximum when 'IGM+' is included. All values are from the literature except the following: The metal content of the ISM at both high 
and low-z, and the baryon content of the ISM at high-z, and the stellar metal content at high-z have been estimated as described in the 
text. The metallicity of the postulated group/field gas at low-z has been estimated to be ~ 0.1 Zq but this has not yet been measured. 
The bottom line shows the baryon density inferred from the CMB measurements (Netterfield et al. 2002) , and the metal content of the 
universe at high and low-z predicted from the integrated star formation history (Pettini 1999; Pagel 2002), with a yield of 0.021-0.042. 
Refs: " Pukugita, Hogan & Peebles (1998). Pettini (1999). = Netterfield et al. (2002). Pagel (2002) 



sity we have associated with the ISM in Table 1 is really a 
minimum. 

The results in Table 1 and Fig. 4 show a number of 
different things. First, there was 6 — 45 times more dust 
in galaxies at 2: = 2.5 than there is today (assuming the 
dust properties used to estimate the mgisses are similar at 
low and high-z, most of the uncertainty here comes from 
the integration of the local DMF). Secondly, as illustrated 
in Fig. 4, the amount of metals traced by dust at high-z is 
3-10 times the total amount of metals in DLA and Lyman 
forest systems. The Lyman forest represents the IGM at 
high-z, while the role of the DLAs is less clear. They may 
be the ISM of galaxies (Wolfe et al. 1986), either the outer 
parts of galactic disks or their precursors (Efstathiou 2000; 
Boissier, Peroux & Pettiui 2002), or dwarf galaxies (Yanny 
& York 1992). Some fraction of them may also be linked to 
outflows from star-forming galaxies and therefore could be 
classed as IGM (Shaye 2001). Whichever is the case, at high- 
z the metal content of the ISM is clearly much greater than 
that of the IGM. At low redshift the situation is reversed: 
the IGM now contains far more metals than the ISM of 
galaxies. This is due to two processes: the consumption of 
ISM gas and metals by star formation and the expulsion 
of metal enriched gas from galaxies by winds. Much of the 
IGM today may be contained in groups of galaxies (this is 
shown in Fig. 4 and Table 1 as 'IGM-I-'). The amount of 



this gas is very uncertain but could be enough to provide 

the baryonic matter necessary to match that inferred from 
recent CMB measurements (Netterfield et al. 2002). The 
amount of metals which may be contained in the postulated 
group IGM gas is even more speculative. The hot IGM gas 
in rich clusters has a metallicity of ~ 0.3 Zq (Renzini 1997) 
but the metallicity of gas in less dense environments, such 
as in groups, may be lower (as predicted by Cen & Ostriker 
1999). We have used an estimate of Z = 0.1 Zq for the 
'IGM-f-' metallicity but stress the large uncertainty in both 
the amount of gas and metals in this category. 

We have also included in Table 1 an admittedly un- 
certain estimate of the total amount of metals contained in 
stars at ^; ~ 2.5. This was estimated from the following three 
populations: 

(i) Lyman Break Galaxies (LBG) - we took the z = 3 
LF from Shapley et al. (2001), corrected to this cosmology 
{nm = 1), and integrated to 0.2L*. Using (B - V) = 0.16 
and M/Lb = 0.15 (Pettini et al. 2001) we used the inte- 
gral of the LF to estimate the stellar mass density p* = 4 x 
lO'^ M©Mpc-^ h75. We used a metallicity of Z = 0.1-0.32® 
(Pettini et al. 2001) to produce an estimate of the metal 
density in stars in bright LBG: p2,,(LBG) ~ 7.5 x lO'' — 
2.3 x 10^ M0 Mpc"'' h75. An obvious uncertainty here is the 
correction to the luminosity density required for dust extinc- 
tion, which could be a factor 2-5 
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(ii) Submm selected sources (> 3 mjy) - the dust mass 
density (from integrating the dust mass function down to 3 
mJy), and the closed box model were used to estimate the 
mass of metals in stars (and the mass in stars). As above, 
we assumed that the galaxies were halfway through their 
evolution, with a gas fraction of ~ 0.5. This gave a stellar 
mass density of p, = 3.8 x 10^ M© Mpc~'^ hys and a stellar 
metal density of p2,.(sub) ~ 2.4 x lO'^ M© Mpc"^ hrs. The 
main uncertainties here are the temperature of the dust (es- 
timates will increase by a factor 2 if the temperature is 20 
K rather than 25 K), and the broadness of the maximum 
in the ratio of dust mass/total mass (Fig. 5). This means 
there is a considerable range in the gas fraction and hence 
the stellar mass associated with a given dust mass. Using 
a model with inflow or outflow will lead to one inferring a 
greater baryonic mass from a given dust mass 

(iii) The rest of the stellar mass at z — 2.5 predicted 
from the SFH and not in bright LBG or SCUBA sources. 
We assumed that this is in lower luminosity systems and 
DLAs and has Z = 0.07 Zq (the average DLA metallicity 
(Pettini 1999)). 

The sum of these three populations is designed to be equiva- 
lent to the total stellar density predicted from the SFH, and 
has a mean metallicity of ~ 0.2 Zq. Correction of the LBG 
values for dust would mean less stars in category (iii) and 
more stellar metals overall (as LBG have higher metallici- 
ties). A comparison of the metals traced by dust (ISM) and 
those in stars at the two epochs leads to a third conclusion: 
at low redshift the ratio of metals in the ISM to those in 
stars is only < 0.08, while at 2; ~ 2.5 this ratio is ~ 4 — 9. 
The chemical evolution model suggests that the ratio will be 
as high as this when only ~ 20 — 30% of stars in the universe 
have formed, in agreement with the star formation history 
inferred from observations. Taken together, this means that 
most metals at z = 2 — 3 are in the ISM of galaxies. This 
is as expected in the early stages of galaxy evolution; later 
as the gas is consumed, the metals become locked into stars 
and expelled into the IGM. 

Table 1 shows that we can now account for all the met- 
als predicted to exist (from integrating the SFH) at both 
low and high redshift, and thus there is no longer a 'missing 
metals' problem. The DLAs do not seem to trace the bulk of 
the metals within galaxies at high-z, suggesting that there 
are indeed biases present in DLA samples against regions of 
high density and obscuration (Pei & Fall 1995) . A survey of 
DLA's selected in the radio (Ellison et al. 2001) has found 
that there is no evidence for a large missing population of 
high column density absorbers. The data allow for QohA to 
increase by a factor of 2 at most. However, the statistics of 
the survey are still somewhat limited, especially at the high- 
est column densities where one may also expect to find the 
greatest quantities of metals (in an unbiased survey). We 
have also shown that the 'missing metals' can be accounted 
for without resorting to a large reservoir of metals in a hot, 
difluse and as yet undetected part of the IGM at high red- 
shift, as has been suggested by various authors (Pettini 1999; 
Pagel 2002). However, there is still a 'missing baryons' prob- 
lem, since the total amount of directly observed baryons is 
less than the amount inferred from recent microwave back- 
ground (CMB) measurements (Netterfield et al. 2002). This 
shortfall is currently thought to be made up of inter-galactic 




Gas Fraction 



Figure 5. Models of how the dust mass of a galaxy changes 
as it evolves. The y-axis shows the predicted dust mass divided 
by the total baryonic mass of the galaxy. The x-axis shows the 
fraction of the baryonic mass that is in the form of gas; as the 
galaxy evolves, it should move from left to right. Spiral galaxies 
today have a gas fraction of 0.1-0.2, eUiptical galaxies have a gas 
fraction closer to zero, but these simple models break down at 
such low gas fractions and can not give accurate predictions for 
the dust to baryonic mass ratio of elliptical galaxies today. The 
models are a closed-box (solid) in which no gas is either lost or 
accreted by the galaxy, a limiting inflow model in which gas is 
accreted (dotted), and two outflow models, where gas is lost at 
a rate proportional to either 1 or 4 times the star formation rate 
(dashed, dot-dashed). Further details are in Eales & Edmunds 
(1996). The models require assumptions about the fraction of 
metals bound up in dust grains (17 = 0.4) and the global yield 
(p = 0.021). 




Figure 6. Predicted dust mass for a spiral galaxy versus gas 
fraction. In this flgure, the dust mass has been divided by the 
dust mass it has at a gas fraction of 0.1. The curves are the 
predictions of the closed-box and outflow models described in the 
caption of Fig. 5. 



gas at both low and high redshift (FHP) ('IGM-I-' in Table 1 
and Fig. 4). 
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5 DUST MASS EVOLUTION AND THE 
SPHEROID CONNECTION. 

The chemical evolution models can also be used to try to 
explain the evolution in the dust mass function shown in 
Fig. 3. Fig. 6 shows how the ratio of dust mass at any point 
in the evolution of a spiral galaxy compares with that today. 
Present-day spiral galaxies have a gas fraction of ~ 0.1 — 0.2, 
and so in Fig. 6 we have shown how the mass of dust in a 
galaxy is predicted to change relative to the dust mass it 
has at / = 0.1. The figure shows that there is relatively 
little scope for the evolution of the dust masses of spirals; 
they can have at the most dust masses which are factor of 4 
greater in the past than their dust masses today. The local 
DMF is composed entirely of spiral galaxies, and using pure 
luminosity evolution (or, in this case, dust mass evolution) 
requires an increase in the dust masses of the local galax- 
ies of a factor ~ 10 in order to reproduce the high redshift 
objects with Ma > . Due to the different shapes of the 
high and low redshift DMFs, the discrepancy is not as large 
for the lowest mass objects, and so the dust masses of the 
faintest submm sources (from the statistical analysis of op- 
tically selected starbursts in the HDF) are consistent with 
simple evolution from the local DMF. Thus the high redshift 
DMF cannot be reproduced by a chemical evolution model 
applied to local spiral galaxies, which argues that the high 
redshift dust sources found in the deep SCUBA surveys are 
a population which is not represented in the dust mass func- 
tion today, i.e. that the high redshift sources are in fact the 
ancestors of the gas and dust poor elliptical galaxies in the 
local universe. 

Present-day elliptical galaxies have a much lower typi- 
cal gas fraction ~ 0.01 and if Fig. 6 were to be drawn for 
ellipticals (showing how dust mass in the past changes rel- 
ative to that at a gas fraction of 0.01), we would see curves 
of similar shape to those in Fig. 6 but which peak at a mgix- 
imum of between 10 and 35 times more dust in the past. 
However, we note that because the gas fractions and dust 
masses of present-day ellipticals are poorly known, and be- 
cause the ratio in Fig. 6 becomes very sensitive to / at low 
gas fractions^ these simple models are unable to predict ac- 
curately the factor by which the dust mass of ellipticals could 
have been greater at their peak dust content. By the time 
elliptical galaxies have reached their current evolutionary 
status, there will be competition between grain formation 
in low-mass giant stellar atmospheres (~ 0.8 — IMq) and 
the very occasional SNae-la, and long-term erosion of grains 
by hot X-ray gas. These arc very different conditions from 
the grain formation in massive and intermediate-mass stellar 
envelopes, SNae-lI, mantle growth in interstellar clouds and 
ISM gas conditions that were present at the peak of their 
dust content. Until a comprehensive model can be based on 
reliable measurements of the gas and dust content of present- 
day ellipticals, it is probably unwise to try to extrapolate the 
dust masses ellipticals should have possessed at their peak 
dust content (robust value) to the dust masses they should 
have today (very uncertain), other than to accept that there 



§ For example, if the gas fraction today were 0.001 then the max- 
imum increases in dust mass allowed become 50 - 250. 



are good reasons why the dust masses could have been much 
greater in the past (by factors 10-100). 

The highest ratio in Fig. 6 is for an outflow proportional 
to 4 times the star formation rate. This is quite extreme, and 
leads to an original mass 5 times greater than the remnant 
today, which would be difficult to produce in most hier- 
archical models of galaxy formation. Outflow episodes are 
expected in a galaxy's lifetime but at a much lower rate, 
or for a shorter time. Table 1 also shows that the change 
in ISM baryons from high redshift to today is of the same 
order as the increase in stellar baryons. This means there is 
not scope for 70-80% of the ISM to have outflowed between 
z ~ 2.5 and today. Some level of outflow is probably nec- 
essary to produce sufficient evolution in the dust masses of 
ellipticals, and to provide an obvious source of at least some 
of the inter-galactic metals observed in the Lyman-forest at 
high redshift, and the inter-cluster medium today. However, 
the elementary outflow models can give rise to problems in 
interpreting the nucleosynthetic yields and metallicities. 

We now use the results so far to address the hypothesis 
that SCUBA galaxies are the ancestors of today's massive 
(L > L*) spheroids. First, we compare the baryons con- 
tained in the two populations. At low-z, we construct a bj 
optical spheroid luminosity function (LF) from the 2dF red- 
shift survey (Folkes et al. 1999) and type speciflc bulge/disk 
ratios (FHP). We then multiply by the stellar mass-to- light 
ratio (FHP) (M/L), = 4.5 - 8.3 to give a spheroid bary- 
onic mass function (assuming that in today's spheroids, the 
stars dominate the baryonic mass). This is integrated down 
to to give the mass density in L > L* spheroids as 
Psph.b = 8.5 - 15.6 X 10'' M© Mpc"^ hrr,. At high-z, we inte- 
grate the DMF down to a mass equivalent to an 850pim flux 
of 3 mjy (^ M^). We then use the ratio of Mdust /Mbaryon 
from the closed box model (Fig. 5) at a gas fraction of 
0.5 to estimate the baryon content of the submm sources, 
Psub.b ~ 7.6 — 15.3 X lO' Mq Mpc~^ hrr,. Since we use the ra- 
tio at its maximum (which is very broad) we are being con- 
servative in our estimate; using outflow or inflow would also 
increase the inferred Mbaryon- The baryonic content of the 
SCUBA sources at high-z is therefore similar to the baryonic 
content of today's massive spheroids, the distant submm 
sources simply being observed at a time when a greater frac- 
tion of those baryons were in the form of gas rather than 
stars. There is a further interesting result, namely that if we 
compare the baryons at z = 2.5 associated with the SCUBA 
sources to those in DLAs, we get the same ratio (1.2 - 2.4) 
£is when we compare the baryons in L > L* spheroids and 
all disk galaxies today (FHP). This coincidence supports the 
idea that spiral disks are the descendants of the high-z DLA 
population (Wolfe et al. 1986; Pettini 1999), and that the 
spheroids are the descendants of the SCUBA sources. That 
there are fax more metals in the SCUBA sources than the 
DLAs is evidence that the bulk of stars in spheroids formed 
at earlier times than those in disks. 

We can also use the chemical evolution model to predict 
the dust mass function that should have been produced by 
the spheroids when they had formed ~ 50% of their stars. 
We again take the local spheroid LF and multiply by the 
stellar mass/light ratio to give a spheroid mass function. 
Using the closed-box model from Fig. 5, we then multiply 
the spheroid masses by the maximum ratio of dust /baryonic 
mass to give a dust mass function. This is shown as the 
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shaded region in Fig. 3. Including outflow has Uttle effect 
on this shaded region, as the decrease in the maximum 
dust/baryonic mass ratio (Fig. 5) is offset by the increase 
in overall baryonic mass at higher gas fractions. The agree- 
ment between this very simple prediction and the observed 
DMF is rather good, supporting the idea that the descen- 
dants of SCUBA galaxies are not the dusty objects in our 
low-z DMF, but the now dust-poor spheroids. Putting this 
another way, bright submm selected sources (> 3 mjy) can 
account for all the predicted dust from L > L* spheroids. 



6 CONCLUSIONS 

Wo have used submillimotro surveys of the local and distant 
universe to produce dust mass functions. Taken together 
these show the evolution in dust mass and ISM metals from 
redshift zero to a ~ 2.5. In conjunction with a simple chem- 
ical evolution model we have found that: 

(i) There was much more dust in galaxies in the past. 
Simple dust mass evolution in local dusty galaxies to match 
the high-z dust mass function is not feasible as the chemical 
models do not permit spiral galaxies to have had x 10 more 
dust at an earlier point in their history. This argues that the 
local submm sources are not the descendants of the sources 
in the deep surveys. 

(ii) Most metals (> 70%) at a ~ 2.5 are found in the ISM 
of galaxies - preferentially in the most obscured and active 
objects. This solves the 'missing metals' problem, which was 
evident when comparing the metals observed in DLAs and 
the Lyman forest with the metal density inferred from inte- 
grating the star formation history. This supports the sugges- 
tion (Pei & Fall 1995) that DLA samples are not unbiased 
tracers of star forming galaxies at high-z, but instead sample 
the lower density /unobscured (and therefore un-enriched) 
environments. 

(iii) Comparison of baryons in L > L* spheroids today, 
and Md > MJ submm sources at high-z, show good agree- 
ment, indicating that the submm sources arc plausibly the 
ancestors of the spheroids. Using the local optical spheroid 
LF and the chemical evolution model, we can also predict 
the dust mass function the spheroids should have had at 
the time when they had formed ^ 50% of their stars. This 
agrees remarkably well with the observed dust mass func- 
tion at high-z, and again strongly supports the hypothesis 
that SCUBA sources (particularly brighter ones) are ellipti- 
cal galaxies in the process of forming most of their stars. 
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tions on improving the manuscript and also Jon Davies, 
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